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We give a short review of recently obtained results on a new lattice formulation of the non-linear 
change of variables which was once called the Cho-Faddeev-Niemi decomposition in SU(2) 
Yang-Mills theory. Based on this formulation, we proposed a new gauge-invariant definition 
, of the magnetic monopole current which guarantees the magnetic charge quantization. We also 

demonstrated the magnetic monopole dominance in the string tension in SU(2) Yang-Mills theory 
on a lattice. Our formulation enables one to reproduce in the gauge-invariant way remarkable 
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1. Introduction 

It is interesting to study color confinement meciianism in Quantum Chromodynamics (QCD). 
The dual superconductor scenario of the QCD vacuum may be a candidate for that mechanism. 
In paticular, it is known that the string tension calculated from the abelian and monopole parts 
reproduces well the original one, once we perform an abelian projection in Maximally Abelian 
(MA) gauge. It is so-called "abelian and monopole dominance". But it has been diffcult to see 
these phenomena in any other gauges. One of the reason of this gauge problem may come from the 
fact that they used the monopoles on the lattice defined by DeGrand and Toussaint(DT) [^. 

Recently, we have demonstrated that the gauge-invariant magnetic monopole can be con- 
structed in the pure Yang-Mills theory without any fundamental scalar field. The success is achieved 
based on a new viewpoint proposed by three of us [Q] for the non-linear change of variables 
(NLCV), which was called Cho-Faddeev-Niemi (CFN) decomposition |^[|], see also We 
have found that the magnetic charge of our lattice magnetic monopole is perfectly quantized. More- 
over, we have confirmed dominance of our magnetic monopole in the string tension, while it was 
first shown in [Q] in the conventional Maximally Abelian (MA) gauge [^. Therefore we can show 
the gauge invariance of the dual superconductor scenario of the QCD vacuum. 

In this paper, we summarize the recent results on a lattice formulation of Yang-Mills theory 
based on NLCV and a gauge-invariant magnetic monopole on a lattice[|l|]. 

2. Lattice CFN variables or NLCV on a lattice 

We have proposed a natural and useful lattice formulation of the non-Unear change of variables 
(NLCV) in Yang-Mills theory corresponding to the CFN decomposition [^,0]. It is a minimum 
requirement that such a lattice formulation must reproduce the continuum counterparts in the naive 
continuum limit. In this stage, therefore, it is instructive to recall how the CFN variables are 
defined in the continuum formulation. We restrict the following argument to SU(2) gauge group, 
for simplicity. 

In the continuum formulation a color vector field n{x) = (?i^(x)) (A = 1,2,3) is intro- 

duced as a three-dimensional unit vector field. In what follows, we use the boldface to express the 
Lie-algebra i'M(2)-valued field, e.g., n(x) := ?ia(-^)7a> 7a = ^Oa with Pauli matrices Oa (A = 1,2,3). 
Then the 5'M(2)-valued gluon field (gauge potential) A^(x) is decomposed into two parts: 

A^(x)=V^(x) + X^(x), (2.1) 

in such a way that the color vector field n(x) is covariantly constant in the background field V^(x): 

= ^^[V]n(x) := d^n{x)-ig[\^{x),n{x)i (2.2) 

and that the remaining field X^(x) is perpendicular to n(x): 

n{x) -X^ix) = 2tr(n(x)X^(x)) = 0. (2.3) 

Here we have adopted the normalization [^{TaTb) = \5ab- Both n(x) and A^(x) are Hermitian 
fields. This is also the case for V^(x) and X^(x). By solving the defining equation ( |2.2[ ), the V^(x) 
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field is obtained in tlie form: 



V^(;c) = c^{x)n{x) - ig ^ [d^n{x),n{x)]. 



(2.4) 



On a lattice, on the other hand, we introduce the site variable , in addition to the original 
link variable Ux.ji which is related to the gauge potential A^(x) in a naive way: * 



Ux,n =exp{-i£gA^{x)), 



(2.5) 



where e is the lattice spacing and g is the coupling constant. Here is Hermitian, nl = rix, and 
Ux.^i is unitary, ?7j^ = U^^^. The link variable Ux.^ and the site variable nx transform under the 
gauge transformation II [0] as 



Ux,^i D.xUx.^^l+^ = U'^^^, nx ^.xtix^l = n!x- 



(2.6) 



Suppose we have obtained a link variable Vx,^ as a group element of G = SU (2), which is 
related to the 5m (2) -valued background field V^(;ic) through 



Vx,n =exp(-/egV^(.x)), 



(2.7) 



where V^(jc) is to be identified with the continuum variable (lA) and hence Vx,^ must be unitary 



A lattice version of (E^) and 



|) is respectively given by 



and 



(2.8) 
(2.9) 



By solving the defining equation ( |2.8[ ) and ( |2.9| ), the link variable V^.^ is obtained up to an 
overall normalization constant in terms of the site variable n.v and the original link variable Ux^^ = 
exp(— /egA^(;t:)), just as the continuum CFN variable and A^(x) in (2.4)^: 



\^{nxUx,^Vl^)=0. 

Both conditions must be imposed to determine V^,^ for a given set of n,- and Ux^^- 



Vx,n = Vx^^[U,n] = Ux,^+nxUx,^nx+^. 
Finally, the unitary link variable Vx,^ is obtained after the normalization: 



Vx,^ = VxAU,n] := Vx^^/xl -tr[V,^^ V,,^] 



(2.10) 



(2.11) 



It is easy to show that the naive continuum hmit e — > of the link variable (2.1 1) reduces to the 
continuum expression (2.4). 



*In general, the argument of the exponential in ( 2.5 ) is the line integral of a gauge potential along a link from x to 
x + fl. Note also that we define a color vector field n(x) := «/i(x)74 in the continuum, while := n^o^ on the lattice 
for convenience. 



See ||l|] for the detail. 
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Therefore, we can define the gauge-invariant flux, @p[U,n], (plaquette variable) by 

%,^v[t/,n] := £-2arg(tr{(l + n,)t/,-,^V,+/i,vVj+v,M^.!v}/tr(l)). (2.12) 



It is also shown that the naive continuum limit of ( |2.12| ) reduces to the gauge-invariant field 
strength; 

&x.iiv - d^Cv - dyc^ +g^^n ■ {d^n x dyu) = G^v{x), (2.13) 

which plays the similar role that 'tHooft-Polyakov tensor played in describing the magnetic monopole 
in Gergi-Glashow model. 

3. Numerical simulations 

It has been shown that the SU(2) master Yang-Mills theory written in terms of (x) and 
n{x) has the enlarged local gauge symmetry G'^^^^.^^i = SU{2)f^^^^i x [SU {2) /U larger than 
the local gauge symmetry SU {2)f^^^^ := SU (2)|^^^; in the original Yang-Mills theory. In order to 
fix the whole enlarged local gauge symmetry G^^f^f, we must impose sufficient number of gauge 
fixing conditions. 

First of all, we generate the configurations of SU(2) hnk variables }, U^,)! = exp [—igeA^ (x)] , 
using the standard Wilson action based on the heat bath method. Next, we define the new Max- 
imal Abelian gauge (nMAG) on a lattice. By introducing a vector field iijc of a unit length with 
three components, we consider a functional F^magIU ,n',^,@] written in terms of the gauge (link) 
variable Ux.^ and the color (site) variable defined by 

F„MAG[U,n-a,&] = £tr(l - (3.1) 

Here we have introduced the enlarged gauge transformation: ^Ux,^ '■= D.xUx^^^l^^ for the link 
variable Ux,^ and ®nv := Qx^x^^^l for an initial site variable nl"^ where gauge group elements Q.x 
and &x are independent SU(2) matrices on a site x. The former corresponds to the SU (2)'° gauge 
transformation (A^)®(x) of the original potential: (A^)®(x) = n(x)[A^(x) + ig^^ d^jO.'^ (x) = 
A;i(x) +D^[A]ft)(x) + 0(ft)^) for the Q.x = e^sa>(x) ^ while the infinitesimal form of the latter reads 

n^(x) = n(x) +gn(x) x e(x) = ii(x) +gii(x) x 0(x) for the adjoint [SU{2)/U{\yf rotation 0;, = 

g''«e(A)_ 

After imposing the nMAG, the theory still has the local gauge symmetry SU (2)"^^ := SU (2);^^^;, 
since the 'diagonal' gauge transformation (0 = 6 does not change the vaue of the functional 
FnMAciU ,n',^,&]- Therefore, configuration can not be determined at this stage. In order to 
completely fix the gauge and to determine rix, we need to impose another gauge fixing condition 
for fixing SU {2)^^^^^^^. For example, we choose the conventional Lorentz-Landau gauge or Lat- 
tice Landau gauge (LLG) for this purpose. The LLG can be imposed by minimizing the function 
Fllg[U-M: 

FiLG[t/;^l] =Itr(l-"[/,,^) ^ 1/4 / d'x [(A^)«(x)]2 (e^O), (3.2) 
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with respect to the gauge transformation Q.^ for the given link configurations {Ux^n}- 

In what follows, the desired color vector field is constructed from the interpolating gauge 
transformation matrix 0;^ by choosing the initial value nf''^ = 03 and 

:= @,Oi@l = 4 = tv[0A@x03®l]/2 {A = 1,2,3). 



(3.3) 



We generate the configurations of the color vector field {Hx] according to the method explained 
above together with the configurations of SU(2) hnk variables {Ux,^.}- Then we can construct 



{Vx^n \U from ( 2.1 1 ). The numerical simulations are performed on an 8"* lattice at j3 =2.2, 2.3, 
2.35, 2.4, 2.45, 2.5, 2.6 and on 16^ lattice at j3 = 2.4 by thermalizing 3000 sweeps respectively. 

3.1 Quantization of magnetic charges 



We construct the gauge-invariant field strength ( [2. 12[ ) to extract configurations of the magnetic 
monopole current {kx,^ } defined by 



k^{s) 



1 



1 



(3.4) 



This definition agrees with our definition of the monopole in the continuum (divided by 2n). 

In our formulation, we have only the real variable 0p[?7,n] at hand, and we are to calculate 
the monopole current using the final term in (^^). Therefore, it is not so trivial to obtain the 
integer-valued k^{s) from the real-valued &p[U,n]. To check quantization of the magnetic charge, 
we have made a histogram of K{s,Ijl) := 2nk^{s) = ^£^vpadv^pa{x + jJ.), i.e., magnetic charge 
distribution. Note that K{s,iJ.) should become a multiple of 2?! if the magnetic charge is quantized. 
Our numerical results show that K{s,iJ.) is completely separated into or ±2?! within an error of 
lO^'*', see Table We have checked that the data in Table [I] exhaust in total all the configurations 
N = 4xS^ = 16384, because the number Ni of links in the J-dimensional lattice with a side length 
L is given by Ni = dU'. This result clearly shows that the new CFN magnetic charge is quantized as 
expected from the general argument. We have observed that the conservation law of the monopole 
current holds, since the number of +2?! configurations is the same as that of —271 configurations. 

Table 1: Histogram of the magnetic charge (value of K{s,lJ.)) distribution for CFN monopoles on 8^ lattice 
at j3 =2.35. 



Charge 


-7.5 
~-6.5 


-6.5 
~-5.5 


-5.5 
~-4.5 


-4.5 
~-3.5 


-3.5 
~-2.5 


-2.5 
~-1.5 


-1.5 
~-0.5 


-0.5 . 
~0.5 


0.5 
~1.5 


1.5 

~2.5 


2.5 
~3.5 


3.5 
~4.5 


4.5 
~5.5 


5.5 
~6.5 


6.5 
~7.5 


Number 





299 

















15786 

















299 






3.2 Monopole dominance of the string tension 

In order to study the monopole dominance in the string tension, we proceed to estimate the 
magnetic monopole contribution {Wm{C)) to the Wilson loop average, i.e., the expectation value of 
the Wilson loop operator (Wf{C)). We define the magnetic part Wm{C) of the Wilson loop operator 
Wf{C) as the contribution from the monopole current kn{s) to the Wilson loop operator: * 



Wm{C) = exp{27ii'£k^{s)N^,{s)), 



(3.5) 



''"The Wilson loop operator Wf{C) is decomposed into the magnetic part W,„(C) and the electric part We{C), which 
is derived from the non-Abelian Stokes theorem, see Appendix B of [0], In this paper, we do not calculate the electric 
contribution (We(C)) where We{C) is expressed by the electric current = dyF^iv. 
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N^{s) =Y^A-^\s-s')^-e^a{i^aS'^^{^ + fi), d'pSj^{s)=J,{s), (3.6) 

.v' 

where A^^ {s) is defined through the external electric source (s) which is used to calculate the static 
potential: d' denotes the backward lattice derivative d'^f{x) = /(x) — f{x — jx), Sp^{s) denotes a 
surface with the closed loop C on which the electric source J^{s) has its support, and A^^{s — s') 
is the Lattice Coulomb propagator. We obtain the string tension by evaluating the average of 
from the generated configurations of the monopoles {k^{s)}. 

We calculate the respective potential Vi{R) from the respective average (IV,- (C)): 

ViiR) = -log{{WiiR,T))/{Wi{R,T-l))} {i = f,m) (3.7) 

where C = {R,T) denotes the Wilson loop C with side lengths R and T. The obtained numerical 
potential is fitted to a linear term, Coulomb term and a constant term: 

Vi{R) = aiR-ai/R + Ci, (3.8) 

where a is the string tension, a is the Coulomb coefficient, and c is the coefficient of the perimeter 
decay: {Wj{R, T)) ~ exp[— a,/?r — c,(/? + T) + afT /R -\ ]. The numerical simulations are per- 
formed on 8^ lattice at j3 = 2.3,2.4, and 16^ lattice at jS = 2.4, 50 configurations in each case. The 
results are shown in Table ^. In order to obtain the full SU(2) results a/, «/, especially, we used 
the smearing method as noise reduction techniques for 100 configurations on 16^ lattice. 



Table 2: String tension and Coulomb coefficient I 







«/ 






2.3(8^) 


0.158(14) 


0.226(44) 


0.135(13) 


0.009(36) 


2.4(8^) 


0.065(13) 


0.267(33) 


0.040(12) 


0.030(34) 


2.4(16^) 


0.075(9) 


0.23(2) 


0.068(2) 


0.001(5) 



Table 3: String tension and Coulomb coefficient II (reproduced from [|7|) 







«/ 


(^DTm 


(XoTm 


2.4(1 6"^) 


0.072(3) 


0.28(2) 


0.068(2) 


0.01(1) 


2.45(164) 


0.049(1) 


0.29(1) 


0.051(1) 


0.02(1) 


2.5(16"^) 


0.033(2) 


0.29(1) 


0.034(1) 


0.01(1) 



We find that the numerical errors at j8 = 2.3 of 8^ lattice are relatively small in spite of small 
size of the lattice 8"*. Moreover, the monopole part reproduces 85% of the full string tension 
(Jf. The data of 8^ lattice at j8 = 2.4 show large errors. The data of 16^ lattice at jS = 2.4 exhibit 
relatively small errors for the full potential and shows that Om reproduces 91% of Of. In general, 
the monopole part does not include the Coulomb term and hence the potential is obtained to an 
accuracy better than the full potential. 

For comparison, we have shown in Table |3| the data of which has discovered the monopole 
dominance for the first time on 16^ lattice where GoTm reproduces 95% of Of. Here GoTm and 
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aoTm denotes the conventional monopole contribution extracted from the diagonal potential 
using Abelian projection in MAG. In particular, the comparison between the data on 16^ lattice at 
j8 = 2.4 and the data on 16^ lattice at j8 = 2.4 in Table |2| reveals that the monopole contributions 
have the same value between the conventional DT monopole and our magnetic monopole. Thus, 
we have confirmed the monopole dominance in the string tension using our magnetic monopole. 

4. Conclusion and discussion 

In this paper, we have proposed a new formulation of the NLCV of Yang-Mills theory which 
was once called the CFN decomposition. This compact formulation of new variables enables us 
to guarantee the magnetic charge quantization. The monopole dominance has been shown anew in 
the string tension. Thus, the magnetic charge quantization and monopole dominance in the string 
tension are confirmed in the gauge invariant way, whereas they have been so far shown only in a 
special gauge fixing called MA gauge which breaks the color symmetry explicitly. 
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